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A method for fast information retrieval from a probe storage device is considered. It is shown 
that information can be stored and retrieved using the optical diffraction patterns obtained by the 
illumination of a large array of cantilevers by a monochromatic light source. In thermo-mechanical 
probe storage, the information is stored as a sequence of indentations on the polymer medium. To 
retrieve the information, the array of probes is actuated by applying a bending force to the can- 
tilevers. Probes positioned over indentations experience deflection by the depth of the indentation, 
probes over the fiat media remain un-deflected. Thus the array of actuated probes can be viewed as 
an irregular optical grating, which creates a data-dependent diffraction pattern when illuminated by 
laser light. We develop a low complexity modulation scheme, which allows the extraction of infor- 
mation stored in the pattern of indentations on the media from Fourier coefHcients of the intensity 
of the diffraction pattern. We then derive a low-complexity maximum-likelihood sequence detection 
algorithm for retrieving the user information from the Fourier coefficients. The derivation of both 
the modulation and the detection schemes is based on the Fraunhofer formula for data-dependent 
diffraction patterns. The applicability of Fraunhofer diffraction theory to the optical set-up relevant 
for probe storage is established both theoretically and experimentally. We confirm the potential of 
the optical readout technique by demonstrating that the impairment characteristics of probe stor- 
age channels (channel noise, global positioning errors, small indentation depth) do not lead to an 
unacceptable increase in data recovery error rates. We also show that for as long as the Fresnel 
number F < 0.1, the optimal channel detector derived from Fraunhofer diffraction theory does not 
suffer any significant performance degradation. 



I. INTRODUCTION 

Probe based data storage with the information 
stored thermo-mechanically as indentations in a polymer 
medium is demonstrated in [3| . Information is written by 
heating the cantilever tip and applying a force to result 
in an indentation being created in the medium. Informa- 
tion retrieval is achieved thermally: the resistance of an 
integrated resistive heater is monitored as the tip scans 
the surface. This resistance drops naturally with can- 
tilever/medium separation due to the decrease in heat 
transfer from the medium to the heater. In order to 
achieve the data rates required by modern storage appli- 
cations, a probe storage system will typically consist of a 
large array of cantilevers reading and writing information 
in parallel. 

In this paper we consider the retrieval of information 
using the optical diffraction patterns resulting from illu- 
mination of the linear cantilever array with a laser light 
source. The main idea is that different combinations of 
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deflected and un-deflected cantilevers in the array result 
in different diffraction patterns, which can therefore be 
used to store information. Such an approach to read- 
back is potentially interesting since it is fast (the inten- 
sity of light can be captured almost instantaneously with 
a photodiode) compared to the thermo-resistive alterna- 
tive (where we must wait for the thermal sensor to reach 
equilibrium temperature). One of the disadvantages of 
the proposed solution is the increased size of the storage 
device as the optical path to detect cantilever deflection 
will occupy more volume than integrated sensors. There- 
fore the method of information retrieval described in our 
paper targets larger form factors of probe storage, for 
example back-up and archival storage. 

Sensing of cantilever deflection by using optical diffrac- 
tion patterns has been introduced by Manalis et al. 0, Q 
By fabricating a cantilever containing intcrdigitated fin- 
gers a diffraction grating is formed. When the cantilever 
is bent, alternating fingers are deflected while remaining 
fingers arc held fixed. The work is extended to paral- 
lel operation by illuminating multiple interdigital can- 
tilevers and capturing the diffraction pattern from each 
cantilever on a separate detector Others have fabri- 
cated a similar design where the fixed fingers are attached 
to a frame surrounding the cantilever and the authors 
recognized the potential of sensing in wave-vector space 



0. Recently it has been shown that the intensity of the 
principle maximum of a diffraction pattern can be used 
to detect cantilever deflections 0- In this case different 
natural resonance frequencies are used to address each 
cantilever. In Q surface topography has been imaged by 
using a digital holographic microscope to capture ampli- 
tude and phase information of the light reflected from an 
array of cantilevers illuminated with a laser source. 

We present a scheme where the cantilevers themselves 
form an optical grating and the state of the array of can- 
tilevers can be retrieved by analyzing the full diffraction 
pattern, which is captured by an array of detectors. 

In order to retrieve the state of deflection of the can- 
tilevers using low complexity signal processing algorithms 
a special modulation scheme for storing information in- 
side diffraction patterns is developed. Binary user in- 
formation is converted to ternary trits which are then 
mapped to a sequence of indentations to be recorded on 
the medium. We show that by analysing the resulting 
diffraction patterns in the Fourier domain it is possible to 
determine the most-likely sequence of trits in linear time 
using a fast recursive graph-based detection algorithm 
that resembles the Smith- Waterman algorithm used for 
genetic sequence matching Q. 

The paper is organised as follows. In Section |IT] the 
channel model is derived using Fraunhofer diffraction the- 
ory. The principle of information storage inside diffrac- 
tion patterns using the central trits is explained. The 
channel model is then extended to capture the effects 
of electronics noise in the system, the signal-to-noise ra- 
tio is defined which allows us to determine the optimal 
indentation depth. The model is then compared with 
experimental diffraction patterns from an array of can- 
tilevers and a high level of agreement is observed. Two 
detection algorithms are presented in Section lllll a simple 
threshold detector and the optimal maximum-likelihood 
sequence detector and their performance is compared. In 
Section lTVl we consider the effect of the global positioning 
errors that characterise highly-parallel probe storage. It 
is shown that using the techniques introduced in [9| it is 
possible to construct an asymptotically optimal iterative 
sequence detector for a large array of probes. We then 
proceed to study the effects of sub-optimal pit-depth and 
modeling noise on the error-rate of the optical data re- 
trieval algorithms. In Section |V] the main results of the 
paper are summarised. In the appendices a derivation of 
the channel model using Helmholtz-Kirchhoff diffraction 
theory is presented. 



II. CHANNEL MODEL 
A. Storing Information In Diffraction Patterns 

Consider a linear array of N reflective cantilevers with 
sensing probes positioned above an area of the medium 
to which binary information string bo bi ...b^^i has 
been recorded using the thermo-mechanical write pro- 
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FIG. 1. Geometry of the optical readout. 



cess (see Figure [T]) . The value b„ = 1 corresponds to an 
indentation in the medium, 6„ = - to the absence of 
indentation, see [lo| for a review of thermo-mechanical 
probe storage. During the process of reading, individual 
probes are pressed against the polymer surface. As a re- 
sult, the n-th cantilever can be in two states: deflected 
if bn = 1 or not deflected if 6„ = 0. Therefore the set 
of N reflective cantilevers can be viewed as an irregu- 
lar optical grating, the irregularity being determined by 
the data string recorded on the medium. If the reflec- 
tive cantilever array is illuminated by a laser source of 
wavelength A then a data-dependent diffraction pattern 
is created. In the practically relevant case of very long 
cantilevers, the diffraction pattern is two-dimensional - 
the distribution of the reflected electromagnetic field de- 
pends on the diffraction angle 9 and distance R to the 
observation point only, see Figure [TJ 

The fact that different binary strings lead to differ- 
ent diffraction patterns allows one to store information 
directly in the patterns rather than in individual indenta- 
tions. The simplest map from binary data to diffraction 
patterns would be an enumeration of all distinct patterns. 

The optical sensors we used in our experiments can 
only detect the field's intensity, but not the phase. The 
principal fundamental question we need to answer is: how 
many bits of information can be stored in the intensity 
of diffraction patterns resulting from 2^ possible config- 
urations of N cantilevers? 

All the length scales characterising the array are much 
greater than wavelength A, see Section III DI for details of 
the experimental set-up. Therefore, we can use Fraun- 
hofer diffraction theory to calculate the distribution of 
intensity of reflected light, see for a review. 



In the Fraunhofer limit any component of the reflected 
electromagnetic field U can be written as a function of 
angle parameter q = kO, where k = 2-k/X is the wave 
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number: 



N-l 



U{q) = Ciq) J2 eWe- 



-iqnd 



n=0 



where d is the cantilevers' pitch (see Figured]) and 



an) 



^2iksbn 



(1) 



(2) 



is the additional phase gained at the n-th cantilever de- 
flected by the distance 6„s, 6„ £ {0,1} is the n-th bit 
of the information sequence; C{q) is a data-independent 
complex amplitude, 



Ciq) = Ao 



kw^ s'm{qw/2) 
2^ {qw/2) ' 



(3) 



where R is the distance to the observation point, w is 
the cantilever's width, the amplitude Aq is related to the 
intensity of the laser source and reflective properties of 
cantilevers. 

For the sake of completeness, we placed the derivation 
of formula ([IJ along with a careful analysis of applicabil- 
ity conditions in the appendices. 

The intensity I{q) = U{q)U{q) of the reflected light 
can be written as follows: 



m^\c{q)\ 



f{n)e 



-iqnd 



(4) 



with coefficients f{n) given by: 

E^=o'~" 62*'=^*"+''- if n > 



(5) 



fi-n) 



if n < 



where 



tn+p,p = K+p -bpe {-1,0, +1}. 



The variables t„+p^p's can be viewed as digits of a bal- 
anced ternary representation of numbers and will be re- 
ferred to as trits. 

Using Q, it is very easy to check that the following 
two information strings lead to the same intensity pattern 

51 = bobi . . . 6jV-2^Ar-l 

52 = bN-ibN-2 ■ ■ ■ bibo, 

where b = 1 — 6 is the operation of bit inversion. In other 
words the transformation of 

Tn ■ bn tlN-l-n, n = 0, 1, . . . iV - 1 

does not change the intensity pattern. An exhaustive 
numerical check for N = 1, 2, . . . , 10 shows that there is 
no other transformation of the information string which 
leaves the intensity pattern invariant. 

Therefore, we can compute the number of distinct 
diffraction patterns corresponding to 2^ binary informa- 
tion strings as the number of orbits of the transformation 



T/v. If N is odd, all orbits have size 2. Therefore the 
number of distinct diffraction patterns is 2^^^. If N is 
even, the calculation is slightly more difficult as Tn has 
2-/V/2 gxed points (strings which are invariant under Tn) 
as well as orbits of size 2. The final answer is 



■)N~1 



if N is odd, 



ttdistinct patterns 

(N) 



2^-1-^2^/2-1 ifiVisevem 



We conclude that a diffraction pattern created by an TV- 
bit binary strings recorded on the medium contains at 
least Zo(722'-^-i-' = — 1 bits of information. 

How can we retrieve this information? The most 
straightforward way would be to enumerate all patterns 
using {N — l)-bit strings. The optimal receiver would 
then compare the pattern generated by the sensors with 
all 2^-1 model patterns stored in its memory and output 
the index of the model pattern which is the most prob- 
able transmitted pattern given the received one. This 
solution is clearly impractical, as the complexity of the 
optimal receiver grows exponentially with N - the num- 
ber of probes in the array. The rest of the paper is dedi- 
cated to the design of modulation and detection schemes 
which would allow low complexity storage and retrieval 
of information using diffraction patterns. The greatest 
difficulty we must overcome stems from the non-linear 
relationship between the information bits and the corre- 
sponding intensity pattern. 



B. Ternary Modulation Scheme 

The expression (|4]) for the intensity of the diffraction 
pattern is a product of data-independent function |C(g)|2 

and a band-limited function J2n=-N+i fin)^~^'^"''^ ■ 
Therefore, to extract the information-dependent coeffi- 
cients {/(n)} from I{q) we need to sample it at finitely 
many points. Given that there is no resonant relation- 
ship between the cantilevers' width w and the array's 
period d, these sampling points can be chosen not to 
coincide with zeros of |C(g)|2. The sequence of sam- 
ple values I{qi)^ I{q2), ■ ■ ■ can be then normalized by 



l/|C7(gi)|M/|C(g2 



using for example a variable 



gain amplifier. This remark allows us to assume from 
now on that |C(q)|2 = 1 and represent the normalized 
intensity of refiected light as 



JV-l 



Iiq)= J2 /(")^ 



— iqnd 



(6) 



Sampling ^ at 2A^ — 1 uniformly spaced sampling 



points 



2'Km 



(2iV- l)d' 



where m = -(TV - l),-(iV - 2), . . . , -1, 0, 1, . . . , iV - 
2,N — 1, we get a system of 2A^ — 1 linear equations 
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for 2A'' — 1 information-carrying coefficients /(n)'s. This 
system can be solved using a discrete Fourier transform 
(DFT, 

N-l 

^^""^^m^ ^ /fen)e"''"'^ (7) 

m=-N+l 

Note that the sampling scheme chosen above is as eco- 
nomical as possible according to the Nyquist-Shannon 
sampling theorem [l3| : the linear sampling frequency is 
equal to ^^^^^^"^ - just over twice the maximal linear 

frequency contributing to Vmax = ^^2tt^'^ ■ 

We conclude that the problem of extracting the data- 
dependent coefficients {f{n)} from the intensity pattern 
I{q) can be solved at the cost 0{NlogN) using the Dis- 
crete Fast Fourier Transform algorithm (DFFT). 

We are now ready to describe a lossy modulation- 
demodulation scheme which would allow us to perform 
the operations Informations Bit pattern written on stor- 
age medium — !■ I{q) — > Information at a cost which scales 
only polynomially with N . 

We firstly note that since the trit tp_q G { — 1,0,1}, 
*M = q and t^^^ —tp,q for any integer n. This simple 
remark allows us to compute the real and imaginary parts 
of/(n): 

N-l-n 

^{f{n))= E (1 + (cos(2M - J (8) 

p=0 
N-l-n 

(9) 

p=0 

While the real part is a non- linear function of the trits the 
imaginary part is linear and crucially it can be written 
as the sum of the first n central trits 

N-l 

tN-i-p,p ^ bN-i-p - bp, < p < — - — . (10) 
Namely, we have the following elementary observation: 

n-l 

==sin(2fcs)^%-i-p,p (H) 

We propose to store information in central trits since 
they can be easily recovered from the output of the DFFT 
using equation ([TT|) : 

„ ^ ^ (/("+!)) _ ^(/N) (12) 
sin(2fcs) sin(2fcs) 

In what follows we will often use the shorthand nota- 
tion for central trits: tp = <jv-i-p,p- 

The method of storing and retrieving information in 
diffraction patterns is now described with assistance 
from Figure [2j Firstly the binary user information can 
be converted to balanced ternary using a lossless en- 
coder, see e. g. [13] ■ The ternary sequence t = 




FIG. 2. Optical Channel Flow 



{tN~i.o, tN-2,1, ■ • ■ , ijv/2-i,w/2) is then mapped to the bi- 
nary sequence b to be recorded on the medium as follows: 
for n = 0,...,^/2- 1: 

r (0,1) if tjv-i-„,n = -1 

{bN-l^n,bn) = < (0, 0) if tN-l-n,n = (13) 
[ (1,0) if ijV-l-«,„ = +1 

(Note the above map between ternary information string 
and binary indentation sequence leads to the loss of infor- 
mation rate: to store information we use only three out 
of possible four configurations of two cantilevers.) At 
the receiving end, the intensity of the diffraction pattern 
can then be measured and the coefficients /(n) recov- 
ered using DFT. The remaining task is to determine the 
sequence t and convert back to binary to return to the 
user. In the absence of noise the central trits are simply 
given by (|12p . In the presence of noise we must detect 
the sequence of trits - this will be discussed in the next 
section. The final step is to convert the ternary sequence 
back to binary. 

A close examination reveals that the computa- 
tional cost of our modulation scheme is bounded by 
0{Nlog{N)), the most expensive step being the Discrete 
Fourier Transform. The price for the simplicity of our 
scheme is the information rate loss mentioned above: we 
can now only store 3^^^ unique ternary sequences on the 
medium as opposed to 2^ binary sequences - this cor- 
responds to a rate of log2(3)/2 sa 0.79 bits. The loss 
occurs as we ignore the information stored in real parts 
of /(n)'s. What we gain however is the ability to store 
and retrieve information optically using low complexity 
signal processing. 

To modify the described scheme to account for the ef- 
fects of channel noise we need to replace step ([12]) with a 
maximum-likelihood sequence detector. To derive such a 
detector, we need to specify the model of channel noise, 
which is done in the next subsection. 
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C. Noise Model 

The experimental characterisation of the optical read- 
out channel is yet to be carried out. In this paper we 
restrict ourselves to the simplest additive white Gaus- 
sian (AWG) model of noise in the system. We expect 
the important sources of noise to be media imperfection, 
electronics noise and shot noise. It is a matter for future 
investigation to verify whether these can be modelled by 
AWG noise. 

Within the AWG noise model, the received signal is 
modeled by adding a noise term to the output of the 
DFFT: forri = l,...,iV/2: 



where {W^n}^=i is a sequence of independent mean-zero 
Gaussian random variables with unit variance and a is 
the standard deviation of the electronics noise. 

To define signal-to- noise-ratio (SNR), notice that in 
the absence of noise, formula ((T2|) is exact. Therefore, it 
is sensible to define the SNR in terms of the increments 
of the received sequence {i?„} as signals: 



sin(2fcs) tp 

p=0 



(jWn (14) 



R 



n+l 



i?„ = .sin{2ks)tn + a{Wn - VK„^i).(15) 



SNR measured in decibels is: 



SNRdB = lOlogi 



sin(2fcs)^ 
3o-2 



(16) 



It is clear that in order to maximise SNR for a fixed a the 
indentation depth s should be chosen so that sin{2ks) = 
1 - the optimal indentation depth is therefore given by 



^opt 



A/8. 



For a typical wavelength A ^ 600 nm, the above formula 
gives the optimal indentation depth of 75 nm. In reality, 
the indentation depth in probe storage is about 5 nm. It 
may be possible to increase the effective deflection size by 
incorporating some kind of mechanical amplification into 
the cantilever design. In this paper we will simply iden- 
tify the smallest indentation depth for which the optical 
readout is possible without modifying the cantilevers. 



D. Experimental Verification of the Diffraction 
Model 
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FIG. 3. Experimental setup of the optical path. A rectan- 
gular shaped sht and a cylindrical lens create a line shaped 
laser spot on the cantilever array. The array is positioned in 
the back focal plane of the lens. The Fraunhofer diffraction 
pattern is projected on the CCD camera. 
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FIG. 4. Comparison of Q with the experimental diffraction 
pattern created by the array of = 5 undeflected cantilevers. 



Although we ultimately foresee the usage of a one 
dimensional array of very fast photodiodes to record 
the diffraction patterns, here we employ a CCD camera 
(Moticam 2300, Motic, Wetzlar, Germany) with a two 
dimensional array of CCD elements (2048 x 1536 pix- 
els) to ease the alignment procedure. The measurement 
bandwidth is 111 Hz. The frame rate of the camera is, 
however, limited to eight frames per second due to the 
limited data transfer rate to the measurement computer. 
A low pass filter is applied to the measured intensity 
pattern to eliminate spatial noise with frequencies higher 
than the width of the information band. 

The experimental and theoretical diffraction curves are 
shown in Fig. An excellent agreement between the 
theory and the experiment is observed: the relative er- 
ror derived using the L^-distance between the theoretical 
and the experimental curves is 3.0 per cent, which is ap- 
proximately 30.6 decibel. 



To verify the accuracy of the channel model wc com- 
pare the diffraction patterns computed within the Fraun- 
hofer theory against experimental diffraction patterns 
captured using the setup shown in Figure [3] The ar- 
ray consists of A'^ = 5 non-deflected cantilevers of width 
w ~ 13.9 /im whose centres are separated by c? = 20 /im. 
It is illuminated with a laser light source of wavelength 
A = 635 nm and a power of 3 mW. 



III. DETECTION ALGORITHMS 
A. Threshold Detector 

In order to derive the maximum-likelihood threshold 
detector we take the increment of the DFFT outputs as 
in (jI2[) and divide by sin(2A:s) so the received signal ?/„ 
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for n = 0, , 



Vn 



,N/2- 

Rn+l 

sin(2fcs) sin(2fcs) 



1 can be written: 

Rn 



tn + 



V2(T 

sin(2fcs) 



Wn (17) 



Then the task of the threshold detector is to examine 
received signals one-by-one and for each y„ infer the most 
likely recorded trit t„: 

in = argmax Pr(t„ = 

tG{-l,0, + l} 

Applying Bayes' theorem and assuming all trits are 
equiprobable the maximum-likelihood criterion can be 
re- written: 



tn = argmax p(y„|i„ = t) 
te{-i,o,+i} 



(18) 



where p{yn\t) is the Gaussian probability density function 
given by: 



Piyn\tn = t) = (StTCT, 



-1/2 



thd) 



exp 



(19) 



where a thd = V^cr/ sin(2fcs) is the effective standard de- 
viation of noise in formula (|17p . Thus the maximum- 
likelihood criterion can be reduced to the following slicing 
of the signal: 



-1 if y„ < -i 
in = { if-|<z>i< 
4-1 if +1 < Vn 



(20) 



B. Sequence Detector 

The maximum-likelihood sequence detector exam- 
ines the whole sequence of received signals Y = 
(Yi, . . . , where F„ is given by: 



R„ 



sin(2fcs) 



Tn{to 



sm{2ks) 



Wni2l) 



where T„(to, ■ • ■ , in-i) G {^?^, • ■ ■ . +ti} is the sum of the 
first n central trits: 



Tn{to, 



Tl-l 

p=0 



(22) 



The sequence detector then determines the most likely 
sequence T = (Ti, . . . , fN/2) 



T = argmax Pr(T I Y) 

T 



(23) 



The above maximum likelihood problem can be easily 
placed in a classical context: by construction the cen- 
tral trits are independent random variables uniformly dis- 
tributed over the set {—1,0,1}. Therefore, the sequence 
of partial sums Ti, T2, is a random walk on the 



set of integer numbers, [15|. Therefore, ([25]) can be stated 
as follows: find the most likely trajectory of the random 



walker originating at given the set of noisy observations 
of the walkers's position at times 1,2,3,.... 

Once the problem {^^) is solved, the most likely se- 
quence of trits can be easily computed from the sequence 
T since for n = 0, . . . , N/2 - 1 

in — Tn+1 — Tn (24) 

where Tq = 0. Again we can invoke Bayes' Theo- 
rem and assuming all sequences T are equiprobable the 
maximum-likelihood criterion (j23p can be expressed: 

N/2 

T = argmaxp(Y|T) = argmax TT p(r„|T„) (25) 

The conditional probability density function p(y„|T„) is 
given by: 

,-1/2 r {Yn~T)^] 



piYn\Tn = T) = (2nat 



exp 



2(72 

Vlt 



(26) 



Where (Tyu = o j sin(2fcs) is the standard deviation of the 
additive noise appearing in equation (|2ip. Substituting 
into ([25]) we arrive at the following expression for the 
most-likely sequence T: 



iV/2 

= argminV(y„-r„)2 



(27) 



We will now demonstrate that the solution to ([77]) can be 
computed recursively using a graph-based algorithm that 
closely resembles the Smith- Waterman algorithm used 
for genetic sequence matching Q. 

Consider the graph with iV/2 -I- 1 time slices indexed 
by n = 0, ...,A^/2. At time slice n we draw 2?! + 1 
vertices (states) corresponding every possible value of the 
variable T„ (the sum of the first n central trits). Note 
that the time slice zero consists of only the zero state. 
We draw an edge (branch) between the state r„_i and 
the state r„ if r„_i + t = T„ where t e {-1,0, +1} 
is a trit. Thus Tn is connected to at most 3 states on 
the left. If Tn is connected to T„_i on the left we say 
Tn-i G L{Tn). A path of length n -|- 1 is is a sequence of 
states (To, Ti, . . . , T„) such that the states arc connected 
by edges, that is for m = 0, . . . , n — 1, T,„ G L{Tm+i)- 
Thus by construction every sequence T = (Ti, . . . , Tjv/2) 
can be drawn as a unique path of length N/2 + 1 through 
the graph starting from the zero state. A sequence of 
partial sums of the sequence of 3 central trits is drawn 
on an example graph in Figure [5] 

To each state in the graph we assign a metric ^„(T„) 
associated with the probability of passing through state 



Tn 



MTn)^iYn-Tn)^ 



(28) 



Next for every path of length n + 1 through 
the graph (Tq, Ti, . . . , T„) we assign a path metric 
PM„(To,ri,...,T„): 



PM„(ro,Ti,...,T„) 



KTjn) 



(29) 
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FIG. 5. A sequence of partial sums of 3 central trits is drawn 
on the graph 



It can be seen that the maximum likelihood sequence is 
the path through the whole graph of the lowest weight: 

T = arg min PM^/2(0, Ti, . . . , Tjv/a) (30) 

T 

Now we set about constructing the recursion to find this 
maximally likely path. Let Pn{Tn) be the set of all paths 
of length 71 + 1 that end in state T„ then we define the 
surviving path metric at state T„ to be the minimum 
weight of all paths that end at state T„: 

SPM„(r„) = ^mm^ PM„(To, Ti, . . . , T„) (31) 

Now the crucial observation is that all paths ending in r„ 
must have passed through one of the states T„_i e L{Tn) 
so the set Pn{Tn) can be written as the union: 

P„(T„)= U P„_i(T„_i) (32) 

T„_ieL(T„) 

The surviving path metric SPM„(T„) can be expressed 
in terms of the surviving path metrics at the states con- 
nected to T„ on the left: 

SPM„(r„)= min [SPM„_i(r„_i)]+Ai„(r„) 

T„_iGL(T„) 

Thus the surviving path metrics can be computed for the 
whole graph recursively. To find the most likely sequence 
we then traceback through the graph following the most 
likely path. Firstly we find the most likely final state 
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FIG. 6. TER performance of detection schemes 

by choosing the state T7V/2 with the minimum surviving 
path metric: 

TV/2 = arg min SPMAr/2 (7^/2) (33) 

Then we traceback through the graph from the state 
following the path of minimum weight. For n = N/2 — 
1,...,0: 

f„= argmin SPM„(r„) (34) 

T„eL(f,.+i) 

To summarise wc now describe the detection algorithm 
in full: 

Algorithm 1 Maximum-likelihood sequence detector 
Inputs : Yi,...,Y^/2 

Outputs : tN~l,0, • • ■ , tN/2.N/2-l 

1: SPMaiO) = 

2: for 71 = 1 to N/2 do 

3: for Tn = — n to +n do 

4: SPMniTn) = minT„„,gL(T„) (^™n-l(T„_i)) + 

M>i(r„) 
5: end for 
6: end for 

7.- T7V/2 = argmin-r^^^ SPMn/2{Tn/2) 
8: for k = N/2 -I to do 
9: f„ = argmin SPMn{Tn) 

10: tn = Tn+l — T„ 

11: end for 

In Figure [5] the trit error rate (TER) of the threshold 
detector is compared against the TER of the maximum- 
likelihood sequence detector. At low error-rates TER = 
10~^ we see an SNR gain of around 2.hdB. 

The complexity of the sequence detection algorithm 
constructed above is 0{N'^) which should be compared 



8 



with the exponential complexity of the exhaustive enu- 
meration scheme. One final remark on the matter of 
detection is that so far all detection schemes only make 
use of the imaginary part of the Fourier coefhcients. In- 
formation about the central trits is also contained within 
the non- linear real part given by (|8]). It may be possi- 
ble to use the real part of the coefficients as a non-linear 
parity-check code if some algebraic or geometrical struc- 
ture can be identified. 



IV. FURTHER CONSIDERATIONS FOR A 
REAL SYSTEM 

A. Global positioning errors 

Up to this point we have not accounted for the po- 
sitioning errors that are present in any real probe stor- 
age device if the information density is higher than 1 
Tb/in^: to read a series of indentations, the array of 
probes has to be moved to the correct sampling position 
with a nanometer precision, which cannot be done with- 
out incurring positioning errors, see [l6j for more details. 

For an optical system the effect of a global positioning 
error is a global reduction in pit depth s. Let us assume 
that the array of N cantilevers at time i suffers a global 
Gaussian positioning error Ji ~ A/'(0, CTj) where aj is the 
standard deviation of jitter noise. Then the imaginary 
part of the Fourier transform output for n = 1, . . . , N/2 



(i) 



b!>^ = sin(2fcs,;) ^ + crW„, 



(35) 



where as before i?Q — and effective pit depth Si is 
given by: 



Si = cxp I I s. 



(36) 



The pit depth is reduced by a multiplicative factor given 
by the simplest probe storage impulse response [l^ 
p{J) = cxp (— J^/PIU^) where PW is a parameter re- 
lated to pulse width. 

All of the detection algorithms discussed so far take as 
input the DFFT output divided by sin(2fcs), see equa- 
tions pT]) and (PT|) . Now, in the presence of global posi- 
tioning errors, the signal degradation due to the loss of 
pit depth (jitter) is a hidden random variable. It is pos- 
sible however, to estimate the global signal degradation 
using techniques developed in Q. 

Consider the following random variable defined for n = 
0, . . . , N/2 — 1 for a fixed time slice i: 

rl:^ - -RI'Ii - - sm{2ks,)tn + a{Wn+i ~ TU„(j37) 

We will show that signal degradation due to jitter can 
be estimated using the empirical average of the square of 



N 



N/2-1 

E 

ri=0 



(38) 



As it is easy to see from ([57]) . random variables 
are not independent. However, they arc 1- 
dependent meaning that subsequences 



(«)2 „(»)2 



„W2 



and rJl^rW^fc+i 



N/2-1 ^'^^ independent for any p and 
k > \. 1-dependent sequences of random variables (and 
more generally M-dependent sequences for any finite M) 
belong to the class of the so called -0-mixing correlated se- 
quences for which the strong law of large numbers can be 
proved using results of the classical paper [l^ ■ (The cited 
paper established the strong law of large numbers which 
is traditionally proved for sequences of independent ran- 
dom variables to the case of correlated random variables 
provided the correlations between distant members of the 
sequence decay at least exponentially fast.) Namely, it 
can be proved that in the limit — > oo the empirical 
sum (1381) converges almost surely to the corresponding 
expected value: 



RRi 



,(■0 



This mean can be computed explicitly: 



E 



2sin(2fcsj) , 

h ^O" 



(39) 



(40) 



Thus the signal degradation due to jitter can be esti- 
mated by: 



N 



'3 {RR, - 2cr2 



(41) 



Taking A^ to infinity this approximation becomes exact: 



N 



lim Cn = sin(2fcsi) 



(42) 



Thus it is proposed that the detection schemes dis- 
cussed in the previous section can still be used in the 
presence of global positioning errors by replacing sin(2fcs) 
in formulas P7)) and ([?!]) with the estimate of sin(2A:si): 
Cn- In order for this approximation to be accurate we 
need a large number of cantilevers. It may not be prac- 
tical to perform optical readback on hundreds or even 
thousands of cantilevers (since the Fresnel number which 
governs the accuracy of Fraunhofcr approximation grows 
with A^) but it is possible to structure the probe array 
as a large number of rows where each row consists for a 
relatively small number of cantilevers (10 for instance). 
Each row can then be read-back optically and then the 
estimate of signal degradation Cn can be computed using 
the output of the DFFT from all rows. 

In Figure [7] we show results of numerical simulations 
which confirm the effectiveness of the proposed detec- 
tion scheme in the presence of global positioning errors. 
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N=1 0x400, )^650nm, s=U8, a =0.1, w=0.5, 10 blocks simulated 
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FIG. 7. Performance of LLN-based detection schemes 



We have simulated a cantilever array with 400 rows each 
containing 10 cantilevers. One million reads with the ar- 
ray were simulated where at each read the whole array 
suffers a global Gaussian positioning error with standard 
deviation aj = Q.IPW. The two detection schemes dis- 
cussed in the previous section have been evaluated using 
the estimate (|41|) (LLN detectors) and compared against 
the same detection schemes but with perfect knowledge 
of the signal degradation due to jitter sin(2fcsi) ("genie" 
detectors). We see that TER is practically the same in 
both cases. 

We conclude with one final remark on the subject of 
the effect of positioning errors for optical storage. For the 
thermo-mechanical read-back channel the loss of signal 
strength due to weak jitter (ctj << PW) can be exam- 
ined by expanding the impulse response into it's Taylor 
series around J = 0: 



p{J) = cxp 



J2 



PW^ 



OiJ' 



We see that the signal strength decays as . For optical 
read-back the signal decays (assuming optimal pit depth) 
can be found by computing the Taylor scries: 



,i„(^exp(-JV™--))=l-(^).7' 



0(f) 



Thus wc find that for an optical read-back system the 
jitter noise strength is of the order of [J /PWY, which is 
much smaller than [J/ PWY - the strength of jitter noise 
in the thermo-mechanical readout scheme. The physi- 
cal reason for a greater resilience of the optical readout 
scheme to the global jitter is that to the leading order 
positioning errors change the phase of the diffraction pat- 
tern leaving the information-carrying intensity invariant. 



Technology 


Wavelength (nm) 


Optimal pit depth (nm) 


CD 


780 


97.50 


DVD 


650 


81.25 


Blu-ray 


405 


50.63 



TABLE L Optimal pit depths for common laser sources 
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FIG. 8. Performance degradation due to small pit depth 



B. Small pit depth 



In section lll Cl it was shown that the optimal pit depth 
is given by s = A/8. In Table |T] the optimal pit depth is 
given for standard laser sources. The thermo-mechanical 
write process produces indentations of depth less than 
10 nm - far smaller than the optimal pit depth even for 
the latest Blu-ray laser sources. It is thus necessary to 
investigate how the performance of the optimal channel 
degrades as pit-depth becomes smaller than optimal. 

In order to determine how the channel performance 
degrades for small sub-optimal pit depth wc must fix the 
level of noise. For the purpose of this initial investigation 
we assume that if pit depth is optimal the signal-to-noise 
ratio is fairly high (greater than 20dB). This is consistent 
with the pre-high pass filter noise levels of up to 14 dB 
observed in the laboratory [3l although of course in any 
real system it is bound to be lower. In Figure[8]we fix the 
noise strength so that for optimal pit depth the SNR is 
22dB and then plot the TER performance of the sequence 
detector as a function of nanometer-scale pit depth for 
various standard laser source technologies. We observe 
that for a Blu-ray laser source trit error rates of the order 
10^"* can be achieved at pit depths of around 10 nm. 

This initial experiment suggests that it may indeed be 
possible to retrieve information at a low error-rate from 
a probe storage system with small indentations optically. 
The next experiment would be to measure the noise levels 
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FIG. 9. Channel performance degrades as the Fresnel number 
increases. 



for a system with cantilevers deflected by indentations of 
a small depth and determine the achievable trit error 
rate. 



C. The effect of modelling errors 

So far all signal processing algorithms have been de- 
rived working in the Fraunhofer limit, that is, when the 
Fresnel number (see Appendix |D| is much less than one: 

F ~ « 1, 

where V is the distance between the probe array and the 
sensor, see Fig. |Tl]for the illustration. As a concluding 
experiment wc study the degradation of channel perfor- 
mance if the observation point is moved closer and closer 
to the cantilever array so that we gradually move out of 
the region of applicability of the Fraunhofer diffraction 
formula. 

In Figure IHl we demonstrate how the signal process- 
ing algorithms we have derived degrade with the Fresnel 
number. The intensity distribution is computed by nu- 
merically evaluating the Kirchoff diffraction formula (jDl[) 
directly. The resulting pattern is pre-filtered and sent to 
the DFFT signal processor as detailed in the section II. 
Information is extracted from the output of DFFT using 
the detection algorithms derived in the previous section. 
Recall that all these algorithms are based on the assump- 
tion that the Fraunhofer approximation is valid. 

We observe that if F < 0.1 there is virtually no per- 
formance degradation but for F > 0.1 the TER increases 
rapidly: for F = 1 the TER of the sequence detector is 
over two orders of magnitude higher than the TER within 
the Fraunhofer limit. We conclude that in order to use 



the optical data retrieval techniques derived in this paper 
we must ensure the Fresnel number F < 0.1. 



V. CONCLUSIONS 

The paper develops a scheme for storing and retrieving 
information in the thermo-mechanical probe storage de- 
vice using diffraction patterns created by illuminating the 
array of actuated probes. Instead of recovering the infor- 
mation independently by individual probes, we propose 
to retrieve information stored in parallel, by analyzing 
the diffraction patterns created by groups of >> 1 
probes. The scheme we developed is slightly reminis- 
cent of partial response channel equalization used in hard 
drives in the sense that instead of trying to avoid diffrac- 
tion, we detect the most likely information sequence to 
give rise to the observed diffraction pattern. The optical 
readout scheme offers several advantages, over a more 
traditional thermo-mechanical method of reading the in- 
formation off the probe storage device: firstly, optical 
sensing is generally faster than thermo-resistive sensing 
due to long equilibration times for the latter; secondly, 
each bit of information stored in a diffraction pattern is 
effectively spread over a region of storage medium, thus 
making the scheme more resilient to strong instances 
of media or electronics noise influencing an individual 
probe; thirdly, there is evidence that the scheme devel- 
oped in the paper offers a better protection of information 
against global positioning errors - the main contributor 
to performance degradation in the nano-scale probe stor- 
age. 

The principle step which enabled us to construct the 
low complexity optical readout scheme is the design of a 
modulation map between user data and the sequence of 
indentations written on the media. The special property 
of this map is a linear relation between user information 
in the balanced ternary form and the imaginary part of 
the Fourier coefficients of the data-dependence diffrac- 
tion patterns. The key is that the patterns themselves 
are generated by the array of probes actuated over the 
storage medium modulated in accordance with the user 
information. 

In order to investigate the performance of the proposed 
optical readout scheme, a channel model has been derived 
for the diffraction pattern obtained when illuminating 
the cantilever array of a probe storage device and it is 
shown that the relative error between the model and the 
diffraction patterns gathered experimentally is less than 
3%. 

It is then shown that by storing ternary informa- 
tion inside the diffraction patterns it is possible to re- 
trieve this information in the Fourier domain using 
low-complexity detection algorithms. In particular a 
maximum-likelihood detector has been derived that de- 
termines the most likely sequence of trits by finding 
the most likely path recursively on a graph. This de- 
tector out-performs the simple threshold detector by 
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2.5dB. The optimal detection algorithm is a version of a 
two-dimensional Viterbi algorithm reminiscent of Smith- 
Waterman sequence matching algorithm used in compu- 
tational biology. The low-complexity of the derived de- 
tection algorithms can be used to find an optimal trade- 
off between the throughput of the detector and its imple- 
mentation complexity. 

The effects of global positioning errors on an optical 
probe storage system have been studied and it is shown 
that the highly parallel nature of probe storage can be 
harnessed to compute an accurate estimate of the result- 
ing global signal degradation (the law of large numbers). 
The detection algorithms fitted with this estimate per- 
form at the same error-rates as the detectors with prior 
knowledge of positioning errors for a large probe array. 
Furthermore, it is shown that the signal degradation that 
occurs as a result of positioning errors is less severe for 
an optical read-back system than the thermo-mcchanical 
equivalent provided that the amount of probes' deflection 
is optimal and equal to A/ 8. 

Finally, it is found that for levels of noise consistent 
with our initial experiments as well experiments with 
thermo-mechanical probe storage, is is possible to re- 
trieve information at low error rates {TER = 10~^) for 
sub-optimal pit depths of the order s=10 nm. The degra- 
dation in terms of error-rate due to modelling noise has 
been quantified and it has been found that as long as 
the Fresnel number of the system F < 0.1, the effects of 
modelling noise are negligible. 

In the present paper we only deal with the problem of 
storing information in diffraction patterns generated by a 
one dimensional probe (sub-)array. The two-dimensional 
generalization of the optical readout scheme and the com- 
plexity analysis of the associated signal processing algo- 
rithms is an open problem. 

We believe that the scheme constructed in this paper 
can find applications outside the field of data storage. 
After all, the problem we are solving is a particular in- 
stance of a more general question; What can be said 
about the surface which gives rise to a given diffraction 
pattern when sensed by an array of probes of known ge- 
ometry? 

The questions investigated in the paper lie on a bor- 
derline between several fields - data storage, optics and 
information theory to name just a few. To make the pa- 
per as self-contained as possible, we included several ap- 
pendices, the aim of which is to derive the Fraunhofer 
formula for the irregular reflective grid systematically 
starting from the two-dimensional wave equation with 
the appropriate boundary condition;_The derivations are 
based on the material presented in 
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FIG. 10. Geometric set up for deriving KirchhofF formula. 
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Appendix A: Kirchhoff Integral Theorem 

Let U, U' be two solutions to Helmholtz equation in 
domain D bounded by the closed curve S = dD: 

(A + k^)U{r) = = (A + fc2)[/'(r), r e D, (Al) 

where k — ut/c, see Fig. [TO] According to Green's theo- 
rem, 

/ dVdivF^ I dS-F, (A2) 
Jd Js 

where F is an arbitrary smooth vector field, dS = nds, 
where n is the inward unit normal vector to the curve S, 
see Fig. [TUl 

Using Green's theorem and Helmholtz equation. 







/ dV {UAU' -U'AU} 
Jd 



D 



dV {div{UVU') - div{U'VU)} 



I dS-{UVU' -U'VU) 
Js 



We conclude that 



/ dS • (C/VC/' - C/'VC/) = 0, 
Js 



(A3) 



III and m 



for any pair of solutions to (jAip in domain D bounded 
by a closed curve S. Note that the curve S does not have 
to be connected. 



ACKNOWLEDGMENTS 



Appendix B: KirchhofF formula in two dimensions 



This work is partially supported by the European 
Research Gouncil within the FP6 collaborative project 
Probe based Terabit Memory (ProTeM) No. 2005 - 



Let G(r | Tq) be the Green's function of Helmholtz 
equation, 

(A + fc2)G(r |r„)==<5(r-ro), (Bl) 
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where rg is the position vector of point P, see Fig. [TOl 
It solves homogeneous Helmholtz equation in the region 
D \ B^{P), where B^{P) is a ball of radius e centered at 
point P € D. Applying (|X3)) to the region D \ B^{P) 
with U'{r) ~ G{r\ Tj^) we get 

/ dS ■ (C/(r)VG(r | r„) - G(r | r,)VU{r)) (B2) 
Js 

= - I dS- (C/(r)VG(r | r,) - G(r | r,)VUir)) , 

where Se{P) = dB^{P) is a sphere of radius e centered 
at P. We can evaluate the right hand side of the above 
equation in the limit e — > 0: the second term vanishes as 
0(e). The first term is equal to 



- [/(lo) hm / dSVG(r | v^) 
= -C/(ro) hm / dl/AG(r|ro), (B3) 

where we applied Green's theorem to convert surface 
integral to volume integral. Using (|Bip to evaluate 
^^'(r I Tq) and taking the limit e — >■ we find that the 
right hand side of (|B2p is equal to — [/(fq). We conclude 
that 

C/do) = - f dS- {U{r)VG{r \ r^,) - G(r | ro)VC/(r)) 
Js 

(B4) 

In two dimensions, 



Gir\r,)^-4'\k\r-r,\), (B5) 



where 



9 /"^ 

H^^\x)^— e"*(t2 _ l)-i/2rfi (B6) 



is a Hankel function of the first kind. Therefore we arrive 
at the following Kirchhoff formula in two dimensions: 

C/(r„) = {u{v)VH[,'\k\v-v^ |) 

^Hl^\k\v-v^\)^U{v)\ (B7) 



Formula (jB7|) expresses a solution to Helmholtz equation 
in terms of its boundary values. In the context of the 
theory of electromagnetic waves in two dimensions, U (r) 
describes the spatial distribution of any component of 
electromagnetic field strength in a monochromatic wave 
with wavelength 2-k /k. 



Appendix C: Helmholtz-Kirchhoff diffraction theory 
for a reflective strip 




FIG. 11. Geometry of a reflective strip 



point P = [H, 0) contained with the domain D enclosed 
by the closed curve S = Si -\- S2 consisting of the arc of 
a circle of radius Ri around point P (5*1) that intersects 
with the line y = —V and the line connecting these two 
points of intersection (S'2). The boundary Si = M \J A 
is split into M, the reflective strip centered at the point 
(0, ~V) with half-width a = w/2, and A the aperture. 

As i?i — )■ 00 the integral on Si vanishes since the 
Somerfield radiation condition is satisfied (only outgo- 
ing waves fall on Si) and the field distribution at Vj^{P) 
can be expressed: 

U{v,,)='- j dS-{u{v)VH{k\v-v^\) 

A+M 

~H{k\v-v,\)^U{v)] (CI) 



Where dS ~ — is the inward normal on S'2. Let 
[/(*)(r) = Aqc^- - be the field distribution of the incident 
plane wave where k = k{—ey). Let A be the wavelength 
of the light source (so that k = 2tt/X). If the half- width 
of the reflective strip is large compared to the wavelength 
(ka ^ 1), then we can use the Kirchhoff boundary con- 
ditions which state: 

(Kl) The field distribution and its derivative are identi- 
cally zero on A: U\a = and = 

(K2) The field distribution and its derivative on M are 
the same as they would be if M extended over the 
whole of S'2: U\m = ^oe -■- and VC/|a/ = -ikU\M 

If the observation point is far away from the reflective 
strip so that k \ r{x) — r q I'^ ^ then we can use the 
following approximations to H{k | r — r g |) and H'{k \ 



Hik\r-r,\) 



irk I r — Tg 



2»(fe|r-r„|-7r/4) ^q2) 



We now apply Formula (|B7|) to the two-dimensional 
geometry given by Figure [TT] consisting of the observation 



H'{k\v-v^\)^iH{k\v-v^\) (C3) 
Substituting into (|C1|) we arrive at the following integral 
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formula for the field distribution at r, 



■o- 



X 1 



V 

I Lix) - lo I 



2k 



I Lix) - lo I 

^ik\rix)-rg\-iTT/4 



(C4) 



Appendix D: Intensity of light reflected by an array 
of probes in Fraunhofer limit 

We now consider the geometry given in Figure [1] 
consisting of an array of N reflecting cantilevers sep- 
arated by pitch d. The array of cantilevers is posi- 
tioned above a polymer medium to which the binary se- 
quence b = (6o,...,6Ar_i) has been recorded by thermo- 
mechanical probe storage write process [l|. The pres- 
ence of a bit 6„ = 1 causes the n-th cantilever to be 
deflected by the indentation depth s. Thus the centre 
of the n-th cantilever is located at position (i/„,— y„) 
where i7„ = nd — {N — l)d/2 and Vn = V + bnS for 
n = 0, . . . , iV - 1 . The horizontal shift of {N - l)d/2 
ensures that the cantilevers are positioned symmetrically 
with respect to the y-axis. As before each cantilever has 
half-width a = w/2. Illuminating the cantilever array 
with a plane light source and treating each cantilever as 
a reflective strip the field distribution at the observation 
point Tq = [H, 0) is given by: 



N-l 



n=0 
X ( 1 + 



2k 



V 



I Lix) - lo I 



Jk\r{x)-rg\-iTT/4 



(Dl) 



where C/„ = e''°^" is the incident field at the boundary 
corresponding to the n-th cantilever. This integral is suit- 
able for numerical evaluation or it is possible to compute 



it analytically by approximating |r — fgl as follows: 

kHx kbnsV 



k\r{x) 



kR- 



R 



R 



(D2) 



Where R — \/H^ + is the magnitude of the vector 
pointing from the observation point to the centre of the 
cantilever array. This approximation is only valid in the 
Fraunhofer limit, that is, when the Fresnel number is 
much less than one: 



F = 



k{{N -l)d/2 + af 
V 



< 1 



(D3) 



By substituting (|D2|) the integral in formula (|Dip can be 
computed explicitly resulting in the following equation 
for the distribution of the amplitude of reflected light at 
diffraction angle 9 = t'An^^{H/V): 



ka sin(0) 

exp [ikshn{l 4- cos(6')) - in{kdsin{9))] (D4) 



xE 

ri=0 

where ^i{9) = k{V + R + {N - l)(isin(6')/2) - 7r/4 is 
a phase shift that does not depend on the deflection of 
the cantilevers, li 9 « 1 (small diffraction angles) we 
can further simplify (|D4p . by writing it as a function of 
parameter q = k9: 



U{q) « A{e) 



2ka? ( sin((7a) 



TXR 



qa 



N-l 



E an)e-'^"'' (D5) 



ri=0 



where ^(n) = e " represents the additional phase 
gained due to the presence of a deflected cantilever, A{9) 
is a complex amplitude ^-independent modulus, and 9- 
dependent phase, which however does not depend on the 
state of the cantilever array. As we are only interested in 
the intensity of reflected light, the phase of A{9) can be 
set to zero. 
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